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Abstract 

The electron-phonon interaction in monolayer graphene is investigated using density functional 
perturbation theory. The results indicate that the electron-phonon interaction strength is of com- 
parable magnitude for all four in-plane phonon branches and must be considered simultaneously. 
Moreover, the calculated scattering rates suggest an acoustic phonon contribution that is much 
weaker than previously thought, revealing an important role of optical phonons even at low ener- 
gies. Accordingly it is predicted, in good agreement with a recent measurement, that the intrinsic 
mobility of graphene may be more than an order of magnitude larger than the already high values 
reported in suspended samples. 
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Graphene, a two-dimensional (2D) sheet of carbon atoms in a honeycomb lattice, contin- 
ues to attract much attention due to its unique physical properties. Aside from a substantial 
academic interest resulting from the relativistic-like behavior of charge carriers, this material 
is considered very promising in device applications as it has an extremely high intrinsic mo- 
bility, even at room temperature. Although in realistic conditions (i.e., placed on a substrate) 
the mobility tends to decrease significantly due to the presence of additional scattering mech- 
anisms at the interfaces much effort is currently being devoted to eliminate, or at 
least minimize, these effects which are detrimental to graphene transport characteristics. 
Therefore, it is crucial to develop an accurate knowledge of the electron-phonon scattering 
as it determines the ultimate limit of any electronic device performance. The strength of 
electron-phonon coupling is typically estimated using the deformation potential approxima- 
tion (DPA); it has been applied for graphene by a number of authors BIB . When the 
corresponding deformation potential constant was estimated from the transport measure- 
ments, however, the results revealed a discrepancy that is too large to be ignored 
Moreover, a very recent observation of mobilities in excess of 10 7 cm 2 /V-s at T < 50 K in the 
decoupled graphene Q drastically departs from the conventionally accepted values, raising 
serious questions about the current understanding of the intrinsic transport characteristics 
of graphene. A detailed theoretical analysis of electron-phonon interaction beyond the DPA 
is clearly called for. 

In this work, we apply a first principles approach based on density functional theory 
(DFT) to calculate the electron-phonon coupling strength in graphene. The obtained elec- 
tron scattering rates associated with all phonon modes are analyzed and the intrinsic re- 
sistivity/mobility of monolayer graphene is estimated as a function of temperature. The 
results clearly elucidate the role of different branches (particularly, the significance of opti- 
cal phonons and intervalley scattering via acoustic phonons) as well as limitations of DPA. 
The obtained "effective" deformation potential constants suggest the possibility that the 
intrinsic scattering rates may presently be overestimated, potentially due to the complex 
influence of extrinsic mechanisms, including the substrate. A subsequent Monte Carlo cal- 
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culation predicts, in good agreement with Ref. [8|, that the mobility of graphene could be 
more than an order of magnitude larger than the already high values reported in suspended 
samples jlj when presumably all of the extrinsic scattering sources are eliminated. 

A key component of our theoretical study is the evaluation of the electron-phonon inter- 



2 



action matrix elements from first principles by employing density functional perturbation 
theory (DFPT) within the DFT formalism This technique has the advantage of deal- 
ing with an arbitrary electronic state k and phononic state q on an equal footing. Each 
phonon is treated as a perturbation of the self-consistent potential created by all electrons 
and ions. The calculation of the potential change due to this perturbation gives the value 
of the electron-phonon matrix element: 



where \i, k) is the Bloch electron eigenstate with the wavevector k, band index i, and energy 
Eiy, ^y^scF i s the derivative of the self-consistent Kohn-Sham potential [9] with respect 
to atomic displacement associated with the phonon from the branch v with the wavevector 
q and frequency coV, q , and M is the atomic mass. A similar calculation has previously been 
performed in graphene for q at high symmetry points K and T, where Kohn anomalies occur 



fsee below) 
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111. However the scattering rate calculation requires summation over the 

n 

entire first Brillouin zone (FBZ). Using DFT and DFPT (Quantum Espresso software [12]), 
the electronic band structure and the phonon dispersion are calculated for a single graphene 
sheet in vacuum. The computation is performed on a hexagonal unit cell with the two atom 
basis, lattice constant a of 2.42 A and approximately 5 A of vacuum space, based on a 
norm-conservin g ps eudopotential and the generalized gradient approximation/local density 
approximation [13] . A Monkhorst-Pack grid [ijj] of 36x36x1 is used with no offset for 
integration. To determine the scattering rate of an electron due to the interaction with 
the phonons, the matrix elements fifc+qk are extracted from the self-consistent electron- 
phonon calculation on the Monkhorst-Pack grid of q-points that covers the entire FBZ and 
incorporates the points of high symmetry: T, K, and M. 

Figure [1] shows the resulting matrix elements for an electron at the Dirac point [i.e., 
k = (471 /3a, 0)] as a function of the phonon wavevector q. The contribution of all six phonon 
branches is considered, whereas the electronic states (both initial and final) are limited to the 
lowest conduction band. The matrix elements display a pronounced anisotropy everywhere 
except near the FBZ center (q = 0), where they are in qualitative agreement with the DPA 
(i.e., nearly isotropic) and also comply with the group theory analysis of electron-phonon 

-^o\g ac \ ~ Q an d lim 9 ^ |g op | ~ constant for 
acoustic and optical phonons, respectively. Another important feature is the Kohn anomaly 



— a strong coupling of the electrons with in-plane optical phonon branches in points of high 
symmetry [h]] . This is prominently illustrated by the three sharp peaks at three equivalent 
K points in the TO mode. The other Kohn anomaly (LO branch, q at the V point) is not 
as distinct. Even more crucial is that their magnitude can be directly validated against 
experimental measurements. Comparison with available experiments and another ab initio 
approach (GW) clearly shows good agreement, though DFT is known to underestimate the 
electron exchange correlation energy in the presence of Kohn anomaly. Specifically, the GW 
calculation and the experiments (Raman spectroscopy and inelastic x-ray scattering - see 
Ref. O and references therein) provide values for \g\ 2 that are larger than our results by 
about 15 % and 60 % for q at the V (LO) and K (TO) points, respectively. This discrepancy, 
however, does not lead to a substantial difference in the calculated scattering rates, since the 
correction is appreciable only in a very small area of the FBZ close to the symmetry points. 
Accordingly, the scattering matrix elements obtained by DFPT are expected to provide 
sufficient accuracy to calculate the intrinsic transport properties of monolayer graphene. 

Assuming low electron concentrations in the conduction band, the electron-phonon scat- 
tering rate can be readily obtained by using the Fermi's golden rule. Figure [2] shows the 
results of the calculation at room temperature (T = 300 K) based on the DFPT matrix 
elements described above [Eq. (pQ)]. Although the result is plotted specifically for electrons 
with wavevector k along the K-T direction, its directional dependence is minor in the en- 
ergy range under consideration. The total electron scattering rate (from all six branches) 
for E ~ ksT is of the order of 10 10 s _1 . This is lower than one might expect based on 
earlier estimates {2, 4, f|; a discussion on the possible reasons is given below. The scattering 
rates for the out-of-plane phonons (ZA and ZO) are much smaller (at least three orders 
of magnitude) in comparison with the rest of the branches (in-plane modes). In 2D struc- 
tures with an in-plane reflection symmetry (e.g. graphene), only in-plane phonon modes 
can couple linearly to electrons 151 ] . On the other hand, all four in-plane branches con- 
tribute comparable scattering rates at room temperature despite the differences in phonon 
dispersion. Accordingly, optical phonons can play an important role even at low electron 
energies. Another interesting feature is a change of slope in the emission rates of TA and 
LA phonons [Fig. 121(a)]. This occurs due to the onset of zone-edge phonon emission near the 
K points (denoted as q = K for simplicity) leading to intervalley transfer of electrons. As 
shown, intervalley scattering by TA phonons can be important since the emission threshold 
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is relatively low: faJg^ — 120 meV is substantially smaller than the corresponding TO and 
LO phonon energies. 

The obtained scattering rates can be used to predict the intrinsic transport properties of 
graphene. The temperature dependence of the electrical resistivity is of particular interest 
as it is believed to be relatively insensitive to the impact of environment (e.g., substrate) 
j^], potentially enabling direct comparison with measurements in real (non-ideal) samples. 
Following the relaxation time approximation, an estimate of the intrinsic resistivity can be 
expressed straightforwardly as: 



Pi 



e 2 D(E F )^r tot (E F )^ , (2) 



where D(Ep) is the electronic density of states at the Fermi level E F , r to t{Ep) denotes 
the total scattering rate at Ep, and vf = 10 8 cm/s is the Fermi velocity. Figure [3] shows 
the result of this calculation for n = 10 12 cm" 2 . At relatively low temperatures (T < 
200 K), the electrons mostly scatter quasielastically with TA and LA phonons (intravalley 
scattering), leading to a linear temperature dependance of resistivity. As the temperature 
increases (T > 200 K), the contributions of both optical phonons and intervalley scattering 
by TA/LA modes increase. The exponential slope in this region signifies the sensitivity of 
large-energy phonon occupancy to temperature. 

Within the DPA, the slope of the intrinsic resistivity in the linear region is proportional 
to the square of the acoustic phonon deformation potential constant 4j. Using the DFPT 
based results, we estimate a value of the "effective" deformation potential constant to be 
D ac ~ 4.5 eV for the combined contribution of TA and LA phonons. When compared to 



those extracted from the experiments, this number is 



airly close to the recently reported 



1 .8 cV 0], while much smaller than others (29 eV [l], 18 eV Q]). On the theory side, 
lowever, D ac of similar magnitude (2.6 eV) was also deduced from a valence force model 



171 ]. With the estimates under non-intrinsic conditions (i.e., experiments) consistently larger 
than the intrinsic theoretical prediction, one likely explanation is that the influence of the 
substrate on graphene electron transport may currently be underestimated. This is the 
most apparent interpretation of the wide spread between experimental data (7.8—29 eV) 
indicative of additional factors or mechanisms in play. 

The issue of identifying intrinsic scattering characteristics is further examined by compar- 
ing our calculations with the results in the decoupled graphene layers (allegedly, the purest 



form of graphene), where the electron mobilities in excess of 10 7 cm 2 /V-s were obtained from 
measurements at low temperatures T < 50 K [S\. For this, a full-band Monte Carlo simula- 
tion is performed. The model utilizes the complete phonon spectrum and electron-phonon 
scattering rates as determined by DFT/DFPT, described above. The inset of Fig. [3] shows 
the obtained velocity-field curves, where the low-field mobility is estimated from the slope 
in the linear region. The calculation produces a very high number of approximately 5 x 10 6 
cm 2 /V-s at 50 K, in good agreement with Ref. ^. Accordingly, the results clearly indicate 
the accuracy of the DFPT scattering rates, particularly those by TA and LA phonons, as 
the low-field mobility in this temperature range is dominated by the interactions with acous- 
tic modes (see Fig. [3]). Note that the calculated strength of the electron— optical phonon 
interaction was validated earlier in connection with the Kohn anomaly. A corresponding 
estimation at 300 K predicts the intrinsic mobility approaching 10 6 cm 2 /V-s that is also 
much higher than previously thought. 

For practical application, it would be convenient to approximate the ab initio scattering 
rates by a simple analytical model. Presently, the expression commonly used for the acoustic 



phonon scattering is given as 
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DLTE* , (3) 



while 



~) = i( E * ~ fe"o) (iV q + !)© {E k - hLU ) + (E k + hio ) iV q ] (4) 

is adopted for the optical phonon scattering p. Here, v s denotes the sound velocity, p m 
the mass density, iV q the phonon occupation number, and B(x) the Heaviside step function. 
Although these formulae are valid under limited conditions [for example, Eq. ([3]) takes into 
account only the LA phonons in the long wave approximation], the deformation potential 
constants D ac and Do can be treated as effective quantities to determine the contribution 
of all the relevant modes. A similar discussion was made above in reference to the slope of 
the electrical resistivity. Then, the total rate can be represented as a sum: 

(l/r) tot = (l/r) ac>KK + (l/r) acKK , + (l/r) op . (5) 

Of the combined contribution of TA and LA phonons, the first term {Xl^acKK re P re_ 
sents electron intravalley scattering described by Eq. ([3]), whereas the intervalley transfer 



i^/ T )acKK' i s calculated by Eq. (j3J), as there is no distinction between acoustic and optical 
modes near the zone edge. The last term accounts for the interactions with both TO and 
LO modes. Figure H] shows the results with the following parameters at 300 K: {^-/T~) acKK 

- D ac = 4.5 eV; {l/r) ac KK , - D = 3.5 x 10 8 eV/cm, hw = Kuj t k a = 124 meV; and (l/r) op 

— D = 1 x 10 9 eV/cm, ftwo = fiuj^P = 164.6 meV. Clearly, Eq. §5§ not only serves as 
an excellent analytical approximation but also reveals the interplay among different phonon 
scattering processes at a given temperature and electron energy. This simple model is ex- 
pected to remain valid for T > T BG with the Bloch-Griineisen temperature T BG ~ 50 — 60 




In conclusion, our first principles analysis clearly illustrates that all in-plane phonons play 
an important role in electron-phonon interactions in graphene and must be considered for 
transport studies at room temperature. Moreover, the results suggest that the influence of 
the substrate may be more significant than previously understood, clarifying at least in part 
the discrepancies observed in the strength of electron-acoustic phonon scattering. Under 
ideal conditions, it is predicted that the mobility of graphene could reach as high as 10 6 
cm 2 /V-s at room temperature. 
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FIG. 1: (Color online) Electron-phonon interaction matrix elements 
calculated by DFPT for k at the conduction band minimum (i.e., the Dirac point) as a function 
of phonon wavevector q. High anisotropy beyond the long wavelength approximation is clearly 
visible. The presence of Kohn anomalies reported earlier is revealed as the peaks (light color) at 
three equivalent Dirac points for TO and at the center of the FBZ for LO. 



9 




FIG. 2: (Color online) Phonon (a) emission and (b) absorption scattering rates at T = 300 K as 
a function of electron energy E^, as k changes along the K-T direction in monolayer graphene. 
The contributions of all six branches are shown separately. The out-of-plane phonons (ZA and 
ZO) do not play an important role. At very low electron energies {E < ksT), the scattering with 
absorption of TO phonons is dominant. 
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Temperature (K) 

FIG. 3: (Color online) Intrinsic resistivity as a function of temperature. At T < 200 K, the 
quasielastic intravalley scattering with in-plane acoustic phonons is dominant (dashed line). At 
higher temperatures, the contributions of both optical phonons and intervalley scattering by TA 
and LA modes (dashed-dotted line) lead to an exponential growth of the resistivity. The inset shows 
the low-field mobilities at T = 50 and 300 K, as determined from the slope of the drift-velocity 
versus electric field characteristic obtained by a full band Monte Carlo simulation. 
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FIG. 4: (Color online) Total scattering rate as a function of electron energy at T = 300 K. The 
dotted line shows the fitted curve calculated from Eq. ©. Clearly, the role of optical phonon 
scattering and intervalley transfer must be taken into account in most cases judging from the 
discrepancy between the solid (total) and dashed (TA & LA intravalley) lines. 
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